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ABSTRACT 

Aims. We investigate the star formation occurring in a region well below the Galactic plane towards the optical reflection nebula 
ESO 368-8 (IRAS 07383-3325). We confirm the presence of a small young stellar cluster (or aggregate of tens of YSOs) identified 
earlier, embedded in a molecular cloud located near the lower "edge" of the Galactic disc, and characterise the young stellar popula- 
tion. We report the discovery of a near-infrared nebula, and present a CO map revealing a new dense, dynamic cloud core. 
Methods. We used near-infrared JHK S images obtained with VLT/ISAAC, millimetre CO spectra obtained with the SEST telescope, 
and optical V-band images from the YALO telescope. 

Results. This star formation region displays an optical reflection nebula (ESO 368-8) and a near-infrared nebula located about 46" 
(1.1 pc) from each other. The two nebulae are likely to be coeval and to represent two manifestations of the same single star formation 
episode with about 1 Myr age. The near-IR nebula reveals an embedded, optically and near-IR invisible source whose light scatters off 
a cavity carved by previous stellar jets or molecular outflows and into our line-of-sight. The molecular cloud is fully covered by our 
CO(/=1-0) maps and, traced by this line, extends over a region of ~ 7.8 x 7.8 pc 2 , exhibiting an angular size ~ 5.4' x 5.4' and shape 
(close to circular) similar to spherical (or slightly cometary) globules. Towards the direction of the near-IR nebula, the molecular 
cloud contains a dense core where the molecular gas exhibits large line widths indicative of a very dynamical state, with stirred gas 
and supersonic motions. Our estimates of the mass of the molecular gas in this region range from 600 to 1600 M . The extinction A v 
towards the positions of the optical reflection nebula and of the near-IR nebula was found to be A v a 3-4magandA v a 12-15 mag, 
respectively. 

Key words. Stars: Formation - Stars: pre-main sequence - ISM: Clouds - ISM: Individual (ESO 368-8, IRAS 07383-3325) - 
Infrared: stars - (ISM:) dust, extinction - (ISM:) reflection nebulae 



1. Introduction 

Since early years, optical reflection nebulae have attracted much 
attention by their extended morphology and clearly non-point- 
like appearance ( e.g. |Hubble |fP922"l |Barnard |[T9T3] [1927] ~ 



|den Bergh |1966| l. In the vicinity of reflection nebulae, different 
types of objects have been found, many associated with young 
stars and the process of star formation (e.g. emission-line stars, 
Herbig-Haro objects, jets and outflows, etc.). For a considerable 
period, these nebulae were among the best tracers of recent star 
formation ocurring in a region of the sky. 

Subsequently, many infrared star clusters and stellar groups 
have been found during surveys towards directions of optical 
and radio HII regions. Testi et al. ( 1998 1 found sources around 
Herbig AeBe stars, part of them related to reflection nebulae. 
|Soares & Bica | p002 2003 ) studied four low-mass star clus- 
ters in the CMaRl molecular cloud, three of them related to 
optical reflection nebulae. Dutra et al. ( 2003 i and |Bica et ak 
[ ( 2003| > investigate d more than 300 optical and radio nebulae 
with 2MASS ([Skrutskie "et al. |1997| > and found 346 new embed- 
ded star clusters, stellar groups, and candidates. 



Star formation sites are commonly associated with molec- 
ular material belonging to the parental cloud out of which the 
young stars form. As the young stars evolve, they feed from the 
available gas and dissipate it. Embedded clusters emerge from 
molecular clouds owing to the disruptive action of newly bom 
stars (Matzner & McKee 2000). Jets and outflows carve cavities 
into the surrounding molecular material. Scattering off the walls 
of these cavities allows photons from embedded young stellar 
objects (YSO) to escape and be detected as near-IR nebulosity, 
denouncing the presence of one or more YSOs (e.g. Yun et al. 
[T9931 [Whitney et al. 1 1993| [2003 1. 



* Based on observations collected at the ESO 8.2-m VLT-UT1 Antu 
telescope (programme 66.C-0015A). 



Furthermore, since most gas and dust that exist in the Galaxy 
are confined to the Galactic disc, and are most abundant toward 
the inner Galaxy (Clemens, Sand ers, & Scoville |1988| l, it is of 
particular interest to search for and characterise the star forma- 
tion ocurring near the "edges" of the molecular disc, either along 
the Galacic disc in the far outer Galaxy or at large vertical dis- 
tances from the Galactic disc. 

As part of our study of young embedded clusters in the third 
Galactic quadrant (e.g. |Yun et af~||2007| |2009| l, we report here 
an investigation of a relatively unstudied star formation region 
(SFR), seen towards the IRAS source IRAS 07383-3325. The 
position of the IRAS source coincides with that of an optical 
reflection nebula. The nebula is evident in the Digitized Sky 
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Survey images (especially in that corresponding to the blue O- 
print) and was discovered by Lauberts (1982 ), who named it 
"ESO 368-8" and at that time classified it as an emission nebula. 
It is also listed in the catalogues of HII regions and reflection 
nebulae of |Brand et al. | ( |1986) and of Magakian (2003| (and 
identified under the names "BRAN 61", "Magakian 287", and 
"GN 07.38.4"). 



Optical examination of the Palomar and ESO plates ( Brand 
et al. |1986| l did not find any evidence of obscuration towards 
ESO 368-8. Thus, the molecular nature of this SFR was first 



revealed in the CO survey of [Wouterloot & Brand (19891 in 



their survey of IRAS sources beyond the solar circle, where 
IRAS 07383-3325 is listed as a good detection (no. 1039). We 
refer to this molecular cloud as WB89-1039, first mapped and 
probed by us using millimetre spectra, and presented here in 
Sect. |X2] 

A preliminary search for the young stellar content of an SFR 
can be performed using 2MASS images in the Ks-band. This has 
been done by |Dutra et al. | ( |2003| l, who indentified a candidate 
"open cluster (IROC)" (DBSB 15) located about 22" east of the 
IRAS source with major and minor axes of 1.9' and 1.7'. Here, 
using a combination of deep near-infrared (JHKs -bands) and 
optical (V-band) images, we characterise the stellar content of 
this SFR towards the reflection nebula ESO 368-8/IRAS 07383- 
3325 and molecular cloud WB89 1039. We confirm the young 
stellar nature of the candidate open cluster of |Dutra et al. | (2003 ). 
We found a group of stars, spatially coincident with DBSB 15 
and extending about 1' to the east of the IRAS source, involved 
in infrared nebulosity. We thus report the discovery of a near- 
infrared reflection nebula seen in the near-IR images. 

The Galactic coordinates of this SFR are I = 248.01°, 
b = -5.46° which, together with its estimated distance of 5 kpc 



(see Sect. 3.2 1, places this region well below the Galactic plane 
in the region of the Galactic warp ( Moma ny et al. ||2006| [May] 
|et al. |[T997| ). Thus, this study adds to what has been our rela- 
tively poor knowledge of star formation occurring in "frontier" 
environments, in this case, the lower "edge" of the Galactic disc. 

Section 2 describes the observations and data reduction. In 
Sect. 3, we present and analyse the results. Section 4 discusses 
the star formation scenario in this region. A summary is given in 
Sect. 5. 



2. Observations and data reduction 

2.1. Near-infrared observations 

Near-infrared (/, H and Ks) images were obtained on 2000 
November 11 using the ESO Antu (VLT Unit 1) telescope 
equipped with the short-wavelength arm (Hawaii Rockwell) of 
the ISAAC instrument. The ISAAC camera dMoorwood et al. I 



1998| contains a 1024 x 1024 pixel near-infrared array and was 
used with a plate scale of 0. 147 arcsec/pixel resulting in a field of 
view of 2.5 x 2.5 arcmin 2 on the sky. For each filter, 6 dithered 
sky positions were observed. A series of 15 images with indi- 
vidual on-source integration time of 4 seconds was taken in the 
/-band. Similarly, series of 20 images, each of 3 second integra- 
tion time, were obtained in the H and the Ks bands. 

The images were reduced using a set of own IRAF scripts 
to correct for bad pixels, subtract the sky background, and flat- 
field the images. Dome flats were used to correct for the pixel- 
to-pixel variations of the response. The selected images were 
then aligned, shifted, trimmed, and co-added to produce a final 
mosaic image for each band JHKs ■ Correction for bad pixels 



was made while constructing the final mosaics that cover about 
3.8 x 3.1 arcmin 2 on the sky. 

Point sources were extracted using daofind with a detection 
threshold of 5cr. The images were inspected to look for false de- 
tections that had been included by daofind in the list of detected 
sources. These sources were eliminated from the source list and 
a few additional sources added in by hand. Aperture photometry 
was made with a small aperture (radius = 3 pix, which is about 
the measured FWHM of the point spread function) and aperture 
corrections, found from bright and isolated stars in each image, 
were used to correct for the flux lost in the wings of the PSF 
The error in determining the aperture correction was < 0.05 mag 
in all cases. A total of 900, 860, and 969 sources were found to 
have fluxes in /, H, and Ks, respectively, and errors < 0.15 mag. 

The JHK S zeropoints were determined using faint in- 
frared standard stars (S427-D and S121-E of the LCO/Palomar 
NICMOS Table of Photometric Standards; see |Persson et al.| 
[] |1998| l and checked with those obtained using 2MASS stars 
brighter than K s = 13.5 mag from the 2MASS All-Sky Release 



Point Source Catalogue (Skrutskie et al. 2006 Cutri et al 



2003). Photometry errors range from about 0.08 mag for the 
bright stars to 0. 15 mag for the fainter ones. From the histograms 
of the JHKs magnitudes, we estimate the completeness limit of 
the observations to be roughly 20.0 magnitudes in the /-band, 
and 18.5 magnitudes in the H and Ks bands. 



2.2. CO mm line observations 

Millimetre single-dish observations of IRAS 07383-3325 were 
carried out at the 15 m SEST telescope (ESO, La Silla, Chile) 
during three observational campaigns in 2001 May, 2002 April, 
and 2002 December. Starting from the (0,0) reference posi- 
tion corresponding to the celestial coordinates of the IRAS 
source, maps were obtained in the rotational lines of 12 CO(1-0) 
(115.271 GHz), I3 CO(1-0) (110.201 GHz), and 13 CO(2-l) 
(220.399 GHz), C 18 O(l-0) (109.782 GHz), and C ls O(2-l) 
(219.560 GHz). The SEST half-power beam width (HPBW) was 
46" at 110 GHz and 22" at 230 GHz, so that, since the adopted 
grid spacing was 46", the I3 CO(2-l) and the C ls O(2-l) line 
emissions were strongly undersampled. The 12 CO(1-0) and 
I3 CO(2-l) maps are composed of 8 x 9 pointings, while the 
13 CO(1-0), C 18 O(l-0), and C ls O(2-l) maps are smaller (5x5, 
6x5, and 3x3 pointings, respectively). 

The adopted integration times were between t mt — 60 s and 
fim = 240 s for each pointing. A high-resolution 2000 channel 
acousto-optical spectrometer was used as a back end, with a to- 
tal bandwidth of 86 MHz and a channel width of 43 kHz, and 
was split into two halves to measure both the 115 and 230 GHz 
receivers simultaneously At these frequencies, the channel width 
corresponds to approximately 0. 1 1 km s _1 and 0.055 km s _1 . The 
spectra were taken in frequency switching mode, recommended 
to save observational time when mapping extended sources. The 
antenna temperature was calibrated with the standard chopper 
wheel method. Pointing was checked regularly towards known 
circumstellar SiO masers, and pointing accuracy was estimated 
to be better than 5". 

The data reduction pipeline was composed of the following 
steps: i) folding the frequency-switched spectrum; it) fitting the 
baseline by a polynomial and subtracting it; Hi) co-adding re- 
peated spectra obtained at the same sky position, weighing with 
the integration time and the inverse of the system temperature; 
iv) obtaining the main-beam temperature T MB by dividing the 
antenna temperature 7a by the tjmb factor, equal to 0.7 for the 
lines at 110-115 GHz and 0.5 for the other lines (at 220 GHz); 
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The spectrum baseline RMS noise (in TmbX averaged over 
all map positions, has been found to be 0.13 K for 12 CO(1-0), 
0.09 K for 13 CO(1-0), 0.15 K for 13 CO(2-l), 0.04 K for 
C 18 O(l-0), and 0.11 K for C ls O(2-l). 



2.3. Optical V -band observations 

Our optical V-band observations were conducted during 1999 
December using the lm YALCQ telescope, at the Cerro 
Tololo Inter-American Observatory in Chile, equipped with the 
ANDICAM camera. The ANDICAM CCD camera contained 
2048 x 2048 pixels, but read-out problems meant that only about 
1400 x 1024 pixels were usable. On a plate scale of 0.3 arc- 
sec/pixel, it provided a field of view of 7.0 x 5.1 arcmin 2 on 
the sky. Standard procedures for CCD image reduction were 
applied. Source extraction and aperture photometry were per- 
formed using IRAF packages. Transformation to the standard 
system was made using a fit to stars common to our images and 
to the US Naval Observatory Catalog. The results were checked 
using the spectrophotometric data of Fitzgerald et al. (1976). 
Their stars no. 27 and 28 fall in our V-band image and their 
magnitudes are compatible (within 0. 1 mag) with the results of 
our photometry. 



3. Results and discussion 

3. 1 . Optical nebula and stars 

Figure [T] shows the central region of the V-band image towards 
ESO 368-8/IRAS 07383-3325. At optical wavelengths, this re- 
gion of the sky does not look particularly spectacular which jus- 
tifies the relatively little attention that has been devoted to it. At 
the centre of the image, a group of a few stars is seen towards op- 
tical nebulosity. Basically, this optical nebula appears to clearly 
involve 3 stars with similar magnitudes and possibly one or two 
fainter stars. This can be more clearly seen in Fig.|2]where the 3 
stars are labelled A, B, and C. The position of the IRAS source 
lies among this group of optical stars (but does not coincide with 
any) and well within the optical nebula. We name these stars here 
as "the optical nebula stars". On the other hand, the IRAS ellipse 
error (not drawn in the figure) encompasses more than one star. 
Given the large IRAS beams, it is likely that all of these optical 
nebula stars contribute to the IRAS fluxes. 

We note that, in our images, the optical nebula ESO 368-8 
extends for about 1' diametre (extension measured down to a 
level of 8<r above the background sky brightness; down to 3<x, 
the size increases just slightly, to 70"), much smaller than the 
quoted size of 4' by |Brand et al. | ( |1986| ). 

Table 1 gives our photometry of the three brightest opti- 
cal nebula stars in the optical V-band and in the near-infrared 
JHK S -bands. 

3.2. Molecular gas and distance 

CO emission has been detected toward this region confirming the 
presence of a molecular cloud. In Fig. [3] we present the CO (J=l- 
0), 13 CO (J=l-0), and C ls O (J=l-0) spectra obtained towards 
the position of the map where the CO emission is strongest: the 
(+46, 0) position of the map, or about 46" east of IRAS 07383- 
3325. CO lines are detected with LSR velocities in the range 
of Vlsr - 47 - 57 km s -1 . The 12 CO(1-0) spectrum peaks at 



Fig. 1. Central region of the V-band image towards ESO 368- 
8/IRAS 07383-3325 (corresponding approximately to the region 
of the cluster DBSB 15, proposed by |Dutra etal7| ( |2003| l). Notice 
the central group of stars involved in optical nebulosity. This 
region covers about 2.56 x 2.56 arcmin 2 and is centred on the 
IRAS source. North is up and east to the left. 
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Fig. 2. Isophotes of the V-band image towards ESO 368- 
8/IRAS 07383-3325, with labels (A, B, and C) referring to the 
3 brightest optical nebula stars (see text). The contour levels are 
8, 12, 20, 30, and 200<x above the background sky value. The 
IRAS source is at the (0,0) position. 

V LS r = 51.3 km s" 1 , while those of 13 CO, and C I8 lines peak 
at 52.1 and 51.6 km s _1 , respectively; these values represent the 
velocities of the gas component in this molecular cloud. 

The kinematic heliocentric distance of this cloud can be 
derived from the peak radial velocities of the spectra. This 
can be done a pplying the circular rotation model by Brand 



YALO is the Yale-AURA-Lisbon-Ohio consortium jBailyn et al. | |& Blitz | ( |1993| . Using V L sr = 51.3 km \ a distance d =s 



5.0 kpc is obtained. Given the Galactic longitude, this implies 
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Table 1. Photometry of the optical nebula stars 



ID 


Multiplicity 


R.A. 


Dec 


m v 


m Ks 


(H-K s ) 


(J-K s ) 






(2000) 


(2000) 










A 




07 40 14.7 


-33 32 34 


14.90 


12.46 


0.24 


0.65 


B" 


1 


07 40 15.8 


-33 32 39 


15.27 


13.04 


0.29 


0.95 


2 


07 40 15.8 


-33 32 40 


15.27 


12.36 


0.26 


0.76 


C 




07 40 16.1 


-33 32 31 


15.20 


13.66 


0.16 


0.45 



" This star is double, as revealed by our near-IR images. The V-band 
magnitude quoted here corresponds to the combined flux of both stars, 
unresolved in the optical image 



s 




Fig. 3. J=l-0 spectra of CO isotopes towards the position of 
the map where the CO emission is strongest (about 46" east of 
IRAS 07383-3325). The 12 CO(1-0) line peaks at an LSR veloc- 
ity of about 51.3 km s _1 . The C 18 O(l-0) line has been multiplied 
by a factor 6. 



a Galactocentric distance of 1 1 .4 kpc to the cloud, in the outer 
Galaxy. We estimate a distance uncertainty of up to 20% due to 
uncertainties in the rotation curve and possible streaming mo- 
tions. 

Given the Galactic latitude, the numbers above imply a ver- 
tical distance to the Galactic plane (defined by b — 0°) of about 
500 pc, well below it and within the region of the Galactic warp 
( May et al. |1997| l. The spatial location of this molecular cloud, 



both at a relatively large heliocentric distance in the outer Galaxy 
and far below the Galactic plane, suggests that there are not 
many background stars (behind the cloud), as discussed in sec- 
tion [332] 



3.2.1 . Molecular gas distribution: a new dense cloud core 

The contour maps of the 12 CO(1-0) and 13 CO(1-0) integrated 
intensities / = JTmb dv are shown in Fig. |4] The CO emission 
is resolved and confined to the mapped region. Our CO maps 
presented here represent the totality of the WB89-1039 cloud, 
observed with unprecedented resolution. Only lines character- 
ized by a signal-to-noise ratio > 3 were considered as genuine 



emission. The cloud appears slightly elongated approximately 
in the direction east-west. In the eastern part, there is a peak of 
strong emission in a core-like or nucleus region; in the western 
part of the cloud, the emission is narrower and more diffuse in a 
tail-like fashion. Thus, we identify a new dense cloud core that 
we designate by "PY10 dense core". 

In all maps, the CO emission peaks at the (+46, 0) (about 46" 
east of IRAS 07383-3325), and decreases steeply in all direc- 
tions. The shape of the strong CO emission is thus approximately 
circular and centred on the (+46,0) position. We name this po- 
sition the "WB89-1039 cloud centre", or the "PY10 core cen- 
tre". This cloud position is where the gas is warmest and has the 
highest column density. The maps of C 18 O(l-0) and C ls O(2-l) 
are not shown because detection of these lines is practically un- 
resolved, restricted to one map position only, the PY10 centre 
position. 

Figure [5] shows the contour map of the 12 CO(1-0) line 
widths. The line widths across the cloud vary from 1.0 km s _1 , at 
the positions where the gas is more quiescent with narrow lines 
well fitted by Gaussians, to 3.4 km s _1 at the PY10 core cen- 
tre where a peak is located. At the (0,0) IRAS position, the line 
width is 1.8 km s _1 , close or intermediate to the value of the qui- 
escent gas. Thus, interestingly, the line widths at the PY10 core 
centre experience a large increase when compared with neigh- 
bouring positions, implying unusual dynamical activity, possibly 
with energetic embedded sources stirring up the gas. 

Taken together, these results suggest that the star formation 
taking place in this small molecular cloud has lead to the for- 
mation of a group of young stars, the optical nebula stars, at the 
location of the IRAS source. As a result, the molecular gas at 
this position has been spent and dissipated. However, to the east 
there is much denser molecular gas, peaking at the PY10 core 
centre, about 46" east of the IRAS source, and showing signs 
of dynamical activity and stirring, possibly caused by additional 
optical invisible sources, thus indicating current star formation 
activity, as demonstrated below. 



3.2.2. Molecular cloud mass and optical extinction 

A few different methods can be applied to estimate the cloud 
mass using line emission of molecular tracers, in this case the 
CO isotopes. A first estimate can be made using the empirical 
correlation between the 12 CO(1-0) integrated intensity and the 
H2 column density along the line of sight, N,j, at position (z, j): 



Nij(H 2 ) = X Tmb dv 



(cm" 2 ) 



(1) 



where the constant X has been determined empiricall y. We use 
the value of 1.9 x 10 20 crrT 2 (K km s -1 ) -1 derived by Strong & 



Mattox ( 1996 ). To obtain the mass from the column density, we 
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Fig. 4. Smoothed contours of 12 CO(1-0) (solid line) and of 
13 CO(1-0) (dashed line) integrated intensities, in the velocity 
range between 47 and 57 km s _1 , observed in the WB89-1039 
and IRAS 07383-3325 region, containing a dense core (PY10) 
about 46" east of the IRAS (0,0) position. The two I2 CO(1-0) 
lowest contour levels are at 1 K km s _1 (5cr) and 3 K km s _1 . 
Subsequent contours are in steps of 5 K km s _1 . The three 
13 CO(1-0) lowest contour levels are at 0.3 K km s _1 (2.5<x) and 
0.6 K km s and 1 K km s _1 . Subsequent contours are in steps 
of 1 K km s . The SEST beam at 115 GHz is displayed in the 
lower right corner. The IRAS source is at the (0,0) position. 
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Fig. 5. Smoothed contours of 12 CO(1-0) line widths observed 
in the same region as Fig. [4] The lowest contour is at 0.6 km s _1 . 
Subsequent contours are in steps of 0.6 km s _1 . The IRAS source 
is at the (0,0) position. The PY10 dense core centre, at the 
(+46, 0) position, is also where the line width peaks. 



also used the "typical" value of 2.72 ( |Allen|1973] l for the mean 
molecular mass fi. We obtain a mass of M co = 1.6 x 10 3 M G for 
the mapped area. 

Furthermore, the derived N(H2) value of the column density 
towards any map position can be used to calculate the optical 
extinction Ay (see Bohlin et al. 1 1978 1 towards that position. 



For the (0, 0)=ESO 368-8/IRAS position and for the (+46, 0)= 
PY10 core centre p osition, using the re lation N(Hi)/Ay = 9.4 x 
10 20 cm -2 mag -1 (jFrerking etai||l982|>, we obtain A v (0,0) = 



4.5 mag and Ay(+46,0) = 15 mag, respectively. 

Another method commonly used to estimate the cloud mass 
consists of assuming LTE conditions and using both an opti- 
cally thick line and an optically thin line, in this case 12 C O(1-0) 
and 13 CO(1-0), respectively. In this method (see e.g. Pineda 
|et al.|[2008] for details), the peak main beam temperature of 
the 1 CO(1-0) line is used to derive its excitation temperature, 
which, for example, is found to be 12.5 K at the map refer- 
ence position. Then, assuming that excitation temperatures are 
the same for both lines, the optical depth and the column density 
of 13 CO are calculated for each line of sight. Finally, to calcu- 
late the total mass, it is necessary to adopt a H? / 13 CO abundanc e 
ratios. The "typical" value is 5 x 10 5 quoted by 



The obtained mass amounts to Mi 



CO(1-0) 



Dickman| ( |1978 i 

= 6 x W~Mr, 



which is compatible with what was obtained by means of the 
previous empirical method. The remaining discrepancy can be 
explained taking into account that the number of 13 CO(1-0) 
spectra, contributing to the total map mass estimate, is markedly 
smaller than in the case of I2 CO(l -0), due to smaller spatial cov- 
erage and lower S IN, and considering that the LTE approxima- 
tion method seems to typically underestimate the 13 CO true col- 
umn densities (Pad oan et al.|2 000). Finally, we should note that 
these two lines, 1Z CO(1-0) and 1J CO(1-0), map different parts 
of the cloud because of their different optical depths, and also 
that, evidently, the empirical method based on the 12 CO(1-0) 
emission is quite reliable in a statistical way, but for individ- 
ual clouds it may give inaccurate results. Therefore, given these 
considerations, the two mass estimates derived here can be con- 
sidered in good agreement. 

The column density of 13 CO can also be used to estimate 
the optical extinction Ay towards the cloud by applying a linear 
relation: 



ci N( li CO) + c 2 . 



(2) 



Several estimates of the c\ and c 2 parameters can be found in 
the literature, varying with the observed region and with the 
cal ibratio n technique. Here, we adopt the results of Frerking et 
[al ( |l982| l for the Taurus region, c\ = 7.1 x 10 mag cm~ z 
and C2 = 1.0 mag. The obtained visual extinction estimates are 
Ay - 3.6 mag for the (0, 0) position and 12 mag for the (+46, 0) 
position, similar to the values obtained from the 12 CO(1-0) data. 

3.3. The young stellar population 

Figure|6]presents the ISAAC JHK$ colour composite image ob- 
tained towards ESO 368-8/IRAS 07383-3325. The three optical 
nebula stars are seen as relatively blue objects slightly below the 
centre of the image. However, the most striking aspect of this 
figure is the red nebula eastward of the centre of the image. This 
infrared nebula, very bright in the K$ -band and invisible at op- 
tical wavelengths, has the shape of a fan (or of a flame) and is 
analysed in more detail in the next section ( 3.3.1| l. 

In Fig. [7] (top), we show the relative positions of the molec- 
ular material (traced by the contour lines of the CO integrated 
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3.3.1 



uster 



intensity) and of the stellar content of this region as seen in our 
deep K-band image (greyscale). We notice the good spatial co- 
incidence of the IR nebula and the PY10 dense core. The lower 
panel of Fig.|7]displays the relative positions of all astronomical 
sources mentioned here. 

Figures [6] and [7] also suggest a slight increase in the surface 
density of stars in the region surrounding the optical nebula stars 
and the infrared nebula, when compared with the rest of the im- 
age. In addition, there is clearly a large concentration of "red" 
stars in this region (much brighter in the Kg -band than in the J 
or //-bands). These aspects will be quantified in section 
but, taken together, they suggest a possible small stellar c 
embedded in a region of local enhanced extinction. 

3.3.1. The infrared nebula 

In Fig. [H] we present the isophotes of the JHKs -band images, 
zoomed-in to show the details of the extended nebular emission 
seen about 46" east of the IRAS source. This location coincides 
with the peak of the molecular emission and of the CO line width 
map (the cloud centre or PY10 core centre). Thus, we name this 
IR-nebula "PY10-IRN". 

PY10-IRN is very bright and extended in the K$ -band image, 
becoming fainter and shrinking in size as we progress towards 
the H and the J bluer wavelength images. In each of the near- 
IR images, the nebula does not contain a detected point-source 
because the full width at half maximum (FWHM) at any point 
within the nebula is much greater than the FWHM of the real 
point sources measured across the images. This indicates that 
we are not seeing the PY10-IRN nebula illuminator source and 
that there is no direct line-of-sight to it. Instead, the nebulosities 
seen in the images are composed of light scattered off the walls 
of cavities carved in the cloud by jets or outflows and into our 
line-of-sight. 

During the star formation process, the envelope and the 
molecular cloud core where a young embedded source resides 
is excavated by a stellar jet and molecular outflow, which are 
typically most intense during the earliest evolutionary stages of 
YSOs: Class ([Andre, Ward-Thompson, & Barsony |1993| l, and 



Class I (Adams, Lada, & Shu [1987 1. In the Class stage, the 



source still has to accrete the bulk of its final mass. In the Class I 
stage, the source is accreting matter from a circumstellar opti- 
cally thick disc and is also surrounded by an infalling envelope. 
The presence of these cavities allows radiation from the cen- 
tral embedded object to escape via scattering off the walls of 
the cavities. Depending on the inclination angle of the cavity 
axis of symmetry to the line-of-sight, different shapes result for 
the nebulae that have been found associated to young embed- 
ded sources. This can be seen in models of embedded infrared 



nebulae that have been produced by Lazareff, Pudritz, & Morin 
(fT^.IWhitney et al. | ( fl993) , | Whitney et al. | ( |2003) , and |Starh 
|et al. |(|2006|). The shape of our IR-nebula can be compared with 



these models. It most resembles models no. 13 or 14 of |Whitney] 
|et al. | ( |1993| l (their Fig. 4b), which corresponds to an embed- 
ded source in a cilindrical cavity with moderate optical depth 
and an inclination angle (of the polar axis) to the line-of-sight 
of cos i = 0.6. Similarly, the Stark et al. (2006) models suggest 
that there is a curved or conical outflow cavity viewed at about 
60° inclination. 

A simple astrometric analysis of the location of the peak in- 
tensity of the nebula in each of the JHKs images supports the 
idea that we are seeing scattered light. The positions of the peaks 
in the JHKs images are not coincident. Instead they are offset by 
about 0.12" and lie close to a straight line, with the values of the 
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WBB9-10S9 (CO cloud) 



ESO 368-8 (opt neb), 
IRAS 07383-3325 




DBSB 15 (cluster) 



PY10 (dense core) 
PY10-IRN (IR neb) 



Fig. 7. ( top:) Smoothed contours of I2 CO( 1 -0) integrated inten- 
sity superimposed on the A^-band image. The axes give the co- 
ordinates relative to the IRAS source, (bottom:) The spatial posi- 
tion of all the relevant objects: 1) the WB89-1039 CO cloud de- 
lineated by the lowest contour of our CO map, 2) the PY10 dense 
core delineated by the highest CO contour, 3) the DBSB 15 clus- 
ter, 4) the ESO 368-8 optical nebula coincident with the IRAS 
source, and 5) the new PY10-IRN near-infrared nebula. 



offsets correlating with the wavelength (i.e., the //-band peak is 
offset about 0.12" from the Ks -band peak and the /-band peak 
is offset about 0.12" from the //-band peak approximately in a 
straight line). The direction of this line is remarkably similar to 
the direction of the elongation of the K s -band nebula. This be- 
haviour is indicative of scattered light. 

In order to provide an estimate of the relative brightnesses of 
the nebula in the JHKs -bands, we performed aperture photom- 
etry of the nebula, for a circular aperture centred on the loca- 
tion of the peak emission, and an aperture radius equal to the 
one adopted for the true stars in the images. The results are: 
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Fig. 6. J (blue), H (green), and K$ (red) colour composite image towards ESO 368-8/IRAS 07383-3325 located in the region of the 
three optical nebula stars (south of centre). Notice the bright and very red near-infrared nebula (east of centre). The image covers 
about 3.8 x 3.1 arcmin 2 . North is up and east to the left. 



-/(PYio-iRN) - 16.9, //(pyio-irn) - 15.1, and ^s(pyio-irn) - 13.7. 
Assuming that the coulours of the nebula are the same as the 
colours of its illuminator embedded YSO, we obtain (J - H) = 
1.8, (H — Kg) = 1.4. These colours are also compatible with the 
source being in the Class I stage of low-mass YSOs (Lada & 
|Adams |1992) . 

3.3.2. Reddening and the embedded stellar population 

Figure|9]shows the histogram for the observed (H - K$) colours 
of the sources detected in both the H and the K$ -band images. 
In the range of to about 0.8, the curve is close to symmetric 
and most sources here correspond to field main-sequence stars 
(located in front of the molecular cloud). The spread in (H-Kg) 
values, about the peak value of (H - K s ) = 0.35, stems from 
the range of intrinsic colours of main-sequence stars and to low 
values of variable foreground extinction in the lines of sight of 
each source. The mean value of the (H - K s ) colours within 
this peak (for (H - K$) < 0.8) is 0.33, with a standard devia- 
tion of 0.15. The spatial location of these peak "blue" sources 



is seen in Fig. [10] (upper panel) where they are represented by 
blue open circles. Red filled circles, on the other hand, repre- 
sent sources in the red tail of the (H - Ks) histogram, the "red" 
sources. The blue sources are scattered randomly and uniformly 
across the image, whereas the red sources are concentrated in 
the region of the molecular cloud. Taken together, these results 
strongly indicate that most red sources are objects associated to 
the cloud, either located behind the cloud or being embedded in 
the cloud and possibly containing near-infrared excess emission 
from circumstellar material. We argue here that, due to the loca- 
tion of this molecular cloud well below the Galactic plane and 
the sparsity of stars in the Galactic halo, few sources seen in the 
near-infrared images are background sources. 

For the general case, this colour excess of (H - Ks) = 0.8 
(about equal to the mean value of (H - K s ) for the field sources 
plus 3 standard deviations) would just separate the objects into 
two groups: foreground stars and {embedded + background) 
stars. However, given that there are very few background stars 
at this location below the Galactic plane, this reddening effect 
separates foreground from embedded stars and thus effectively 
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Fig. 8. Isophotes of the JHK$ -band images towards a newly 
discovered near-infrared nebula ("PY10-IRN", coincident with 
the PY10 core centre, located about 1' east of ESO 368- 
8/IRAS 07383-3325). The small ticks on the axes indicate the 
size of 1". North is up and east to the left. Contour levels are at 
5, 10, 20, 60, 180, 600, and 2400cr above the sky background. 



selects YSOs embedded in the cloud. Thus, the large majority of 
the red sources are likely to be YSOs embedded in the molecular 
cloud. 

Blue sources, on the other hand, may be composed of a mix 
of foreground field sources and YSOs in a more evolved evolu- 
tionary stage. These more evolved young stars, if present, could 
be pre-main-sequence objects or even intermediate-mass or mas- 
sive main-sequence stars formed in this cloud, which evolve 



150 - 



100 - 




(H-Ks) 

Fig. 9. Histogram of the observed (H - K s ) colours. The well- 
defined peak is composed of foreground field sources ("blue" 
sources with colours (H - Ks) < 0.8). The red wing of the dis- 
tribution is composed of sources spatially concentrated in the 



region where the molecular cloud is present (see Fig. 10 1. 



much faster than their lower-mass siblings formed at the same 
time. Their higher masses would also contribute to their being 
bluer and thus not being told apart by red colours. 

3.3.3. Young stellar objects 

Using the 648 point sources detected in all three J, H, and 
K s -bands, we plotted the near-infrared colour-colour diagram, 
(/ - H) versus (H - Ks ), shown in Fig[TT] Most stars are located 
within the reddening band where stars appear if they are main- 
sequence stars reddened according to the interstellar extinction 
law (Rieke & Lebofsky 1985), which defines the reddening vec- 
tor (traced here for Ay = 10). Pre-main-sequence YSOs, or mas- 
sive main-sequence stars recently formed in this region, which 
have had time to clear the inner regions of their circumstel- 
lar discs, lie in this region as well. Giant stars appear slightly 
above this band. Stars located to the right of the reddening band 
are likely to be embedded young star objects with infrared ex- 
cess emission from circumstellar material ( Ada msTLada, & Shu| 
[p987l l. 

For the sources that lie inside the reddening band, the high- 
est value of (H - Ks) is about 1.14. Using the mean value of 
(H - Ks) = 0.33 for field stars (according to Fig. |9j, we obtain 
a colour excess E(H - K s ) = 0.81 due to intra-cloud extinction, 
which corresponds to about a maximum visual extinction pro- 
duced by the cloud core of about Ay — 13 (Rieke & Lebofsky 
1985). This value is very similar to those derived towards the 
densest part of the cloud (position (+46", 0) or PY10 core) using 
CO spectra (see section [3~.2.2| ). 

The location of the vertical dashed line, derived from Fig. [9] 
is at (H - K s ) = 0.80. The two groups of sources, blue sources 
with (H-K s )< 0.80 and red sources with (H - Ks) > 0.80, are 
very differently distributed on the colour-colour diagram. A large 
fraction of the red sources are located outside and to the right of 
the reddening band, whereas the blue sources mostly occupy the 
inside of the reddening band. Thus, most red sources are likely 
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Fig. 10. ( top:) Spatial distribution of all the sources seen both 
in the H and in the K s band images. Blue open circles represent 
sources with values of (H - K) lower than 0.8. Red filled cir- 
cles represent sources with values of (H - K) higher than 0.8. 
This plot is centred on the IRAS point source, (bottom:) Spatial 
distribution of all the sources seen simultaneously in the /, H 
and /fj band images. Red filled circles represent sources to the 
right of the reddening band (see next figure). Blue open circles 
represent sources located inside the reddening band. 



to be YSOs. Given their spatial concentration (Fig. \TQ[ , these 
red sources, with (H - K s ) > 0.80, seem to represent a small 
young embedded stellar cluster of about 46 young stars form- 
ing in the molecular cloud. The actual number of stars formed in 
this cloud is likely to be larger for at least three reasons. Firstly, 
we chose a conservative value of (H - Ks) = 0.80 (3<r above 
the mean value of the (H - Ks) of field main-sequence stars). 
Secondly, the presence of the optical nebula stars seems to in- 
dicate that there are some young stars in a more advanced stage 
of the star formation process, already free of circumstellar ma- 
terial and exhibiting blue colours, thus not pinpointed by our 
colour selection criterion. Thirdly, there are IR-excess sources 
with (H - K s ) < 0.80 that are YSOs belonging to the cluster. 
The presence of this YSO population is likewise supported by 



Fig. 10 (lower panel) where a plot similar to the upper panel is 



shown. We notice that unlike the upper panel, in the lower panel 
red filled circles represent sources with IR excess emission (de- 
rived from Fig. 111. The distribution of the IR excess stars is 




Fig. 11. Near-infrared colour-colour diagram of the region to- 
wards IRAS 07383-3325. The solid line represents the loci of 
unreddened main-sequence stars (Bessel & Brett 1988}, while 
long-dashed lines indicate the reddening band. The reddening 
vector indicates the direction of the shift produced by extinc- 
tion by dust with standard properties. The location of the vertical 
dashed line, derived from Figj9] is at (H - Ks) = 0.8. The cross 
in the lower right corner represents a typical error bar. 



are concentrated at the location of DBSB 15, between ES0368- 
8/IRAS 07383-3325 and the PY10-IRN nebula (closer to the lat- 
ter), and trace the young stellar population forming in this cloud 
core. Thus, as mentioned above, the specific location of this SFR 
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in the Galaxy makes the two approaches illustrated in Figs 
almost equivalent in pinpointing young stars. 

Given the uncertainty in the number of stars forming in the 
cloud, we can make a rough estimate of the star formation effi- 
ciency of this molecular cloud. Assuming that there are 46 solar 
mass stars and 3 optical nebula stars that could reach about 6 so- 
lar masses each (see next section), we obtain a total stellar mass 
of ~ 64M Q . For a total mass of the gas in the cloud between 
600 and 1600 M , we obtain a star formation efficiency in the 
range of 4-10%. Even though the uncertainty is large, these val- 
ues are similar to the lower end values of star formation efficien- 
cies found in cluster environments within the local star formation 
regions (e.g., L1630; |Lada |1999 1. 



4. Discussion 

The two non-spatially coincident nebulae, ESO 368-8 (seen in 
the optical) and PY10-IRN (seen in the near-IR images), to- 
gether with a small enhancement of the surface stellar density are 
the main signs of star formation in the form of a small embed- 
ded young cluster in the WB89-1039 molecular cloud. Several 
aspects make this region particularly interesting. The first aspect 
is the Galactic location, in the outer Galaxy and well below the 
Galactic plane. Even though not in the far outer Galaxy (defined 
as Galactocentric distances >13.5 kpc; this number is based on 
the radial distribution of CO emission from the studies of Digel 



et al. 



1996 Heyer et al. [1998), this location shares at least 
one aspect with the far outer Galaxy: it is close to the "edge" of 
similar to the distribution of the large (H - Ks) sources. Both the Galactic molecular disc (in this case, the lower edge), about 
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500 pc (vertical distance) below the plane of the Galaxy defined 
by b — 0°. Interestingly, star formation occurs frequently and 
with no major differences in these edge environments, both in 
the far outer Galaxy (e.g. |Brand & Wouterloot ||2007[ |Yasui &| 
Kobayashi |[2008| |Yun et al. ||2009| l, and also well below the 
Galactic plane as shown here. The relatively small number of 
molecular clouds at this large vertical distance ( |May et aL~[l 997 ) 
do not appear to significantly affect the occurrence and proper- 
ties of the star formation process. The only significant difference 
of these "edge" star formation sites is that they do not seem to 
form massive stars (e.g. |Yun et al. 2007, 2009 ). In fact, based on 
the observed fluxes, no massive star is expected to be forming in 
this region. 

In the following, we discuss the nature of the optical neb- 
ula stars and the presence of the two non-spatially coincident 
nebulae, focusing especially on what they suggest for the star 
formation scenario in this molecular cloud. 



4. 1 . The optical nebula stars 

The nature of the stars involved in the optical nebula can be 
investigated using the (H - K S ),K S colour-magnitude diagram 
seen in Fig. 12 In this figure, we have marked the positions of 
the optical nebula stars (A, B, and C). The optical nebula star B 
is in fact a double star resolved in our near-IR images (whose 
stars are named here Bl and B2). The best fit to these sources is 



achieved using the 1 Myr isochrone of Siess , Dufour & Forestini| 
[ ( 2000| l with extinction Ay = 3. This value is in good agree- 
ment with the value of Ay derived towards these stars (the (0,0) 
position) using the CO data (Sect. |3.2.2| i. Bearing in mind the 
relatively large uncertainties when fitting isochrones to colour- 
magnitude diagrams, we estimate that the optical nebula stars 
have an age close to 1 Myr and are intermediate-mass stars (from 
spectral type B5 to A9). Furthermore, using Ay = 3, the absolute 
magnitudes in the V-band of the optical nebula stars span values 
that are a few magnitudes greater than the corresponding val- 
ues for intermediate-mass main-sequence stars. This is compati- 
ble with these stars being pre-main-sequence intermediate-mass 
stars whose luminosities are still higher than the corresponding 
ones when they reach the main-sequence. 



Figure 13 shows the colour-magnitude diagram for the 
sources seen in both the H and the Ks -bands but excluding those 
that exhibit signs of excess emission due to circumstellar ma- 



terial (as inferred from Fig. 111. In this way, the location of a 



star on this diagram depends on its distance (vertical shift), on 
its extinction (shift along the reddening vector) and on its age 
(shift given by pre-main-sequence theoretical models of stellar 
photospheres with no circumstellar emission). The two groups 
of sources (filled red sources with (H - K s ) > 0.8 and open 
blue sources with (H - Ks) < 0.8) are separated by the vertical 
dotted line. The blue sources represent mostly foreground main- 
sequence stars, at variable distances and small extinction. The 
red sources can be best fitted by a 1 Myr isochrone with vari- 
able values of extinction, as is known to exist toward embedded 
sources, between about Ay = 8 and Ay = 15, consistent with the 
values of Ay towards these sources (mostly located around the 
(+46,0) position) derived from the molecular data (Sect. [3~Z2| . 
Thus, this plot supports the idea that we are in the presence of 
a young stellar cluster or aggregate with about 1 Myr age and 
composed of at least some tens of stars. 
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Fig. 12. Colour-magnitude diagram with the positions of the 
optical nebula stars (A, Bl, B2, and C) marked by filled dots. 
The dashed line is the zero-age main-sequence (ZAMS; |Siess, 
|Dufour & Forestini |2000) > for the cloud distance and with zero 
extinction. The solid fine is the 1 Myr isochrone with extinction 
Ay = 3. The reddening vector indicates the direction of the shift 



produced by extinction by dust with standard properties ( Rieke 
|& Lebofsky (1985) . 

4.2. The two reflection nebulae 

The two reflection nebulae, the optical nebula ESO 368-8 and the 
near-IR nebula PY10-IRN, are separated by about 46" (1.1 pc). 
Not surprisingly the location of the optical nebula is relatively 
poor in gas, when compared with the location of the near-IR 
nebula where the amount of gas (and the extinction) is much 
greater. Do the two nebulae mark two star formation sites in the 
same cloud, or can the two nebulae be considered coeval and 
two manifestations of the same episode of star formation taking 
place in the cloud? In the latter case, the two different stages of 
the young stars, the optical nebula stars and the red sources close 
to the near-IR nebula, could stem from different stellar masses. 
The optical nebula stars should be sufficiently more massive than 
the stars forming close to the near-IR nebula. Because of their 
higher masses, they would have evolved faster and dissipated 
the cloud material while their lower mass siblings to the east are 
still heavily embedded and extincted at optical wavelengths. 

On the other hand, if the two nebulae represent two non- 
coeval star formation sites, could this be the case for triggered 
sequential star formation? Using the gas temperatures from the 
CO observations to derive the sound speed in the cloud, we con- 
clude that it would take about 2.3 Myr for a pressure wave to 
travel from the optical nebula to the near-IR nebula (for a pro- 
jected distance of 1.1 pc between the two nebulae). If the optical 
nebula stars have triggered star formation at the near-IR loca- 
tion, they had to have done it at least 2.3 Myr ago. This value is 
of the order or greater than the age derived for the optical neb- 
ula stars. This fact, and the ages derived in the previous section, 
favour the scenario of the nebulae being considered coeval and 
representing two manifestations of the same single star forma- 
tion episode. Moreover, the higher velocities inferred from the 
line widths within the PY10 core, peak at the position of the 
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Fig. 13. Colour-magnitude diagram for the sources seen in both 
the H and the Ks -bands with no signs of excess emission due 
to circumstellar material. The dashed line is the zero-age main- 
sequence for an arbitrary distance. The two solid lines are the 
1 Myr isochrones with extinctions Ay — 8 and Ay = 15, respec- 
tively. Filled circles represent the red sources with (H - K s ) > 
0.8. Open circles represent the blue sources with (H—Ks) < 0.8. 

near-IR nebula. This indicates that the source of gas stirring here 
is likely to be an embedded source at this position, not a source 
located far out at the position of the optical nebula where line 
widths are not enhanced. 

5. Summary 

We present a study of the molecular gas and young stellar popu- 
lation seen towards the region of the reflection nebula ESO 368-8 
and the IRAS source IRAS 07383-3325, in the molecular cloud 
WB89-1039 located in the outer Galaxy and well below the 
Galactic plane. We confirm the presence of a young stellar clus- 
ter (or aggregate of tens of YSOs), identified by |Dutra et al. | 
(2003), composed of low and intermediate -mass stars, with no 
evidence of high-mass stars. 

This star formation region is associated with the ESO 368-8 
optical reflection nebula and with a near-infrared nebula (PY10- 
IRN) located about 46" (1.1 pc) from each other. The two nebu- 
lae seem to be coeval and to represent two manifestations of the 
same single star formation episode with about 1 Myr age. The 
near-IR nebula appears to be composed of scattered light off a 
cavity carved by previous stellar jets or molecular outflows from 
an embedded, optically, and near-IR invisible source. 

The molecular cloud was fully covered by our CO(/=1-0) 
maps and, traced by this line, extends over a region of ~ 7.8 x 
7.8 pc 2 , exhibiting an angular size ~ 5.4' x 5.4' and shape (close 
to circular) similar to spherical (or slightly cometary) globules. 
Towards the direction of the near-IR nebula, the molecular cloud 
contains a dense core (PY10) where the molecular gas exhibits 
the large line widths indicative of a very dynamical state, with 
stirred gas and supersonic motions that could generate turbu- 
lence. 

Our estimates of the mass of the molecular gas in this region 
range from 600 to 1600 solar masses. In addition, the extinction 



Ay towards the positions of ESO 368-8 (optical reflection neb- 
ula) and of PY10-IRN (near-IR nebula), derived by means of dif- 
ferent column density methods, was found to be Ay ^3-4 mag 
and Ay =s 12-15 mag, respectively. We estimate a value of 
~ 4 — 10% for the star formation efficiency of this molecular 
region. 
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